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Abstract 

We study a model of carbon nanotube with a half filled conduction band. At this filling the 
system is the Mott insulator. The Coulomb interaction is assumed to be unscreened. It is shown 
that this allows to develop the adiabatic approximation which leads to considerable simplifications 
in calculations of the excitation spectrum. We give a detailed analysis of the spectrum and the 
phase diagram at half filling, and discuss effects of small doping. In the latter case several phases 
develop strong superconducting fluctuations corresponding to various types of pairing. 

PACS numbers: PACS No: 71.10.Pm, 72.80.Sk 
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I. INTRODUCTION 



Carbon nanotubes have attracted an enormous amount of attention and generated an 
immense level of theoretical and experimental activity. It has been rightly pointed out 
that single-wall carbon nanotubes (SWCN) represent ideal one-dimensional systems and 
one may expect to observe here exotic phenomena characteristic for strong correlations in 
one dimension. Most of the attention has been concentrated on a possibility of Tomonaga- 
Luttinger liquid. Though such liquid is certainly a very interesting object, it is just one of 
many wonders strong correlations can produce. 

In this paper we concentrate on the physics of an armchair SWCN at half filling. The 
features which dominate this physics are the unscreened Coulomb interaction and Umklapp 
processes. The importance of the unscreened Coulomb interaction for nanotubes away from 
half filling was already noticed in the previous studies 

HQ 

. At half filling some additional 
factors come into play. It was realized in that the long-range Coulomb interaction makes 
the operators responsible for the Umklapp scattering terms relevant. Consequently, the gaps 
for collective excitations generated by such operators are not exponentially small, as it would 
be away from half filling, but have power-law dependence on the Umklapp scattering matrix 
elements. Such enhancement of the gaps increases chances for their experimental observa- 
tion. Unfortunately, in their further analysis the authors of resorted to Renormalization 
Group (RG) equations, which for systems with many fields do not provide much insight 
into the properties at the strong-coupling fixed point. RG equations also do not take into 
account the drastic difference between velocities of the plasma modes and all other collective 
excitations. Instead of adding difficulties, however, the difference in velocities leads to con- 
siderable simplifications in actual calculations (see ^). In this paper we shall exploit this 
feature to our advantage and develop an approach based on the adiabatic approximation 
similar to the one used for the problem of electron-phonon interaction. As a result we will be 
able to provide a rather detailed information about the spectrum and the phase diagram of 
the system. As will be shown below, the interplay of the long range Coulomb interaction and 
Umklapp processes at half filling gives rise to a phase diagram which includes several inter- 
esting strongly correlated states. The system has a hidden Z2XZ2XZ2 symmetry and these 
phases are conveniently classified as different symmetry breaking patterns of this group. 



2 



II. THE PROBLEM OF SINGLE-CHANNEL WIRE AT HALF FILLING 



To warm up, let us first recall the problem of the unscreened Coulomb interaction for 
a single chain (or a single channel quantum wire) at half-filling. The charge and spin 
sectors decouple; the Tomonaga-Luttinger liquid Hamiltonian for the charge sector should 
be supplemented by the Umklapp term which contains only the charge field $c : 

V = Ua {R+L^R + l L i + h.c.) = cos[V87r$c], (1) 

where Uao is the 2kp = n Fourier component of the interaction. The full Hamiltonian 
density for the charge sector is TC = 7i + V, where 
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*c) 2 + (9.4c) / dyV(x - y)(A,*c) 

2 

V(x) = v 2 S(x) + (2) 
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v being the Fermi velocity. This model is very close to the sine-Gordon one. The spectrum 
contains solitons and their bound states (breathers or excitons). To get their spectrum one 
can just expand the cosine around its minimum and obtain: 



lu 2 = (vqY 



2e 2 

1 + — m(l/g(V 
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where the breather gap is m& = (2/n)\/UeF, ep = irylao. The breathers are effectively 



optical phonons of the one-dimensional Wigner crystal 
estimated in |4(, a better estimate is 



The soliton gap is larger; it was 



M s = -[1 + (2e 2 /7rv) \n(v/M s a )} JuT F (4) 

7T v 



The spectrum in the spin sector is gapless, and the spin velocity is approximately v. Thus 
already for a simple one-chain problem, the long-distant Coulomb force brings three impor- 
tant new features: (i) strong upward renormalization of the charge velocity seen in Eq.(|3*]l. 
and (ii) the presence of exciton modes in the charge sector, (iii) power law dependence of 
the mass scales M s ,rrib on the Umklapp matrix element U. 

III. THE EXCITATION SPECTRUM OF NANOTUBE AT HALF FILLING 

SWCNs are manufactured by wrapping two-dimensional graphite sheets into cylinders 
(compactification). The electronic spectrum of an infinite graphite sheet contains two Dirac 
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cones with different chirality, centered at different points of the Brillouin zone (the blue dots 
on Fig. 1). Under compactification the spectrum is divided into subbands corresponding to 
different quantized values of the transverse momentum. In the so-called armchair nanotubes 
k y = remains an eigevalue and the spectrum in the lowest subband stays gapless; at half 
filling it is also doubly degenerate at low energies. This degeneracy is a vestige of the double- 
cone structure of the two-dimensional dispersion. We discuss the spectrum in some detail 
in Appendix A; more detailed information can be obtained from the book Q|. 
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FIG. 1: The Bravais lattice cell (shown in dotted lines) including two inequivalent carbon ions 
(shown in red and green) and the Brillouin zone for graphite, ei and e2 are basis vectors of the 
diatomic unit cell. The positions of the tips of two Dirac cones labeled by 1 and 2 are represented 
by blue and brown dots. 

In Appendix B we derive the bosonized form of the effective low-energy Hamiltonian. 
The total Hamiltonian density is 

H = tto[4 +) ] + ~ E [(^ (a) ) 2 + v\d x ^f] + V (5) 

where 7io is given by Eq.(J2J); $c + ^ * s the symmetric charge mode. The label a takes three 
values, a = (c, — ), (s, +), (s, — ), corresponding to the antisymmetric charge field and sym- 



metric and antisymmetric spin fields, respectively (see Appendix A). The interaction density 
V contains Umklapp terms 



V = -or^cos^Trtf 



(+)^ 



x [# c cos(v / 4^$ c ] ) + (93 ~ 9x) cos(V^$i +) ) + 9s cos(V47r9( } ) - g x cos(v / 4^$i } )] 
+••• (6) 

where the dots stand for all other terms that do not involve $c an d stay marginal even in the 
presence of the unscreened Coulomb intreraction. The couplings g c , g\ and #3 are determined 
by the lattice Hamiltonian. In Appendix B we estimate them for realistic carbon nanotubes, 
where, as it turns out, g c = g\ and, hence, there are only two independent coupling constants. 
The unscreened Coulomb interaction strongly reduces the scaling dimension of the operator 
cos(v^4tt ( £ > c + ' ) ) i n the long-wavelength limit, making it smaller than 1. Thus the Umklapp 
terms (jUJ) become strongly relevant, with the scaling dimension almost equal to 1. This 
circumstance dramatically increases the values of the gaps. 

One may well expect that the double degeneracy of the electron band in carbon nanotubes 
makes the problem similar to the problem of two interacting chains (the 'ladder' problem) 
much discussed in literature. This is indeed the case. To put the problem in a broader 
context, let us discuss two extreme types of ladders: (i) the ones where two chains are 
placed far apart such that there is no direct tunneling between them, and (ii) those in which 
the interchain tunneling is stronger than all the interactions. Hamiltonian (J5I6|) is formally 
equivalent to the first case describing two well separated chains. As a matter of fact, the 
Hamiltonian for case (ii) is not very different. The calculations done in [7[ for the ladder 
with strong interchain tunneling yields the same Hamiltonian as ©, but with $c,- field 
being substituted by its dual counterpart, 0c,-: 

V B = V A [$c,- -> ©c,-] (7) 

Therefore one expects that the two models have the same excitation spectrum though the 
response functions are different due to the different field identification. 

Hamiltonian @ is not very convenient to analyse in its bosonic form since the effective 



potential contains mutually nonlocal and noncommuting fields, $ S) _ and O s ,-. The p 
becomes significantly more transparent when one uses the refermionization procedure 
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Fermionizing all the fields except of we obtain 

H = tt [4 +) ] + H-R)d x R f + L\d x Lf) + ^ j2(-X a R d x x a R + xld xX l) 

3 



— cos(V^^ c +) )[ 9c (R\L f - h.c.) + 2g t ]T X rXI + ^sXrxI], (8) 

a=l 

9c = gi, 9t = (g^ - 9i), 9s = -gi - #3, 



where \Ri Xl are Majorana (real) fermions and Rf, Lf are Dirac fermions emerging from the 
fermionization of the operator cos(v^47r$c,-) or, in the second case, cos(\/47r@c,-). Notice 
that Eq.(JHJ) is manifestly SU(2)-symmetric. The triplet of Majorana fermions (a = 1,2,3) 
transform according to the spin S=l representation of the SU(2) group, whereas the fermion 
labeled by a = is a singlet under the SU(2). 

Following Q], we shall handle Hamiltonian (jHJ) using the adiabatic approximation, whose 
validity is guaranteed by the fact that the velocity in the symmetric charge sector is strongly 
enhanced by the long-distant Coulomb interaction with respect to the bare Fermi velocity. 
The results have certain similarity with the SU(4) theory proposed in but also contain 
important differences, which we shall discuss. 

Thus, from the ^E^-mode point of view, the other degrees of freedom are static. Inte- 
grating over this mode, one obtains an effective potential for the fermions in the form of 
the ground state energy of the sine-Gordon model. According to the energy density is 
proportional to the square of the breather mass fib of such sine-Gordon model and in the 
regime of small sine-Gordon coupling constant [3 2 is equal to 

2 3 

E « « -v/{-Kav c )[ig c {R)L f - h.c.) + 2ig t £ X rXI + 2i</.xM> (9) 

c a=l 

where v c ~ Vp[l + (2e 2 /irvp) In^ep/M)} 1 / 2 and v ~ Vp. So as we see, the integration over the 
fast mode gives rise to the mass terms for all of the fermions. The fermionic modes acquire 
gaps: 

M C ,. = ^L, M S ,« = ^L M S ,„ = ^ (10) 
ixav c ixav c irav c 

They correspond to neutral excitations with the quantum numbers given in the table that 
follows; see also Fig.2. As far as the fast modes are concerned, one has to treat differently 
the ones with and without electric charge. The neutral modes do not involve solitons of 
field $ c +) . One can get a good estimate of their spectrum replacing the term in the square 



brackets in expression for the Hamiltonian (JSj) by a constant and expanding around the 
minimum of the cosine potential. As a result one gets the same spectrum as for the single 
chain (J3J) with 

m 2 b = vr(M c 2 „ + 3M 2 4 + M S 2 S ) ln(v/a) 

where with the logarithmic accuracy M is either Mc,- or Ms- Thus 

m b /M ~ [\n(v/aM)] l/2 > 1 

which further supports the adiabatic approximation. However, as we have already men- 
tioned, its validity is already assured by the difference in the velocities. Apart from the 
breather modes there are massive modes corresponding to half-period solitons in <&c + ^ with 
zero modes of the fermions bound to them. These excitations have the largest gap (we call 
it M e ) and carry quantum numbers of electron. 

Table 1. 



Mass Q S V 

M s>t triplet 1 

M s>s singlet 

M c ,_ vortex ±2 

mj breathers 

M e electron ±1 1/2 ±1 

The results for the spectrum are summarized in Table 1.; see also Fig. 2. Excitations 
are characterized by quantum numbers associated with the full continuous symmetry group 
of the effective model ((H). These are the total charge Q measured in units of the electron 
charge e, total spin S, and "vorticity" V. Let us comment on the latter. While the global 
U(l) phase invariance and the spin SU(2) symmetry, leading to conservation of the total 
charge and spin, 




are exact symmetries of the original, microscopic Hamiltonian, an extra ("flavor") U(l) sym- 
metry generated by global phase transformations of the /-spinor (Rf,Lf) (or equivalently, 
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a.) 



b.) 



FIG. 2: A schematic picture of the excitation spectrum: a) at half filling and b) at small doping. 

by uniform translations of the dual field 0£ - ^), emerges in the low-energy limit only. This 
symmetry leads to conservation of the "flavor" charge 



which, together with conservation of Q, implies independent conservation of the particle 
numbers at each Dirac point, — (Q ± V)/2. 



FIG. 3: A staggered flux (orbital antiferromagnet) state corresponding to a nonzero value of 
vorticity V. The arrows indicate local currents flowing across the links of the a (red) and b (green) 
sublattices. 

The nomenclature "vorticity" we have chosen follows from the microscopic origin of the 
flavor density. Indeed, consider a lattice operator describing a current flowing around the 




(12) 




S 



elementary plaquette of sublattice a or b: 



JpiaqMr) = 1 
v = a,b 



^(r)^(r - e 2 ) + if>l{r - e 2 )^(r + ei) + ^(r + ei)Vv(r) - h.c. (13) 



From this construction it follows that the state with a nonzero (j p iaq,a(b) ) would represent an 
orbital antiferromanet, or a staggered flux phase, realized on the sublattice a (b), as shown 
in FigJSJ Projecting (JT3j) onto the k y = subband and passing to the continuum limit, one 
makes sure that the sum j p i aqja + j p i aq ,fo does transform to the flavor charge density 



Pi 



[■ RLR*t ■ - ■ RtaR2a + {R ~+ L)] ~ d x $ 



Previous attempts to study the problem of two coupled channels have mostly relied on 
the assumption of equal velocities with the subsequent use of RG analysis. The approach 
was pioneered by Lin et al. [lfj who argued that at strong coupling the spectrum of the 
two-chain problem at half filling acquires a higher symmetry, such as SO (6) or even SO (8). 
In fact, Gell-Mann-Low equations alone are insufficient to extract information about strong 
coupling regime (they have to be supplemented by Callan-Symanzik equations for the physi- 
cal observables) and therefore cannot provide a legitimate ground for such conclusions. Our 
approach is based on a different assumption; here the long-range Coulomb interaction legit- 
imates a clear separation of scales between single particle excitations and collective modes. 
This significantly simplifies the calculations providing one with the well-controlled approxi- 
mation. It is instructive to compare the results with those conjectured in [llj]. Though the 
structure of the multiplets is the same as in the SU(4)~SO(6) theory, there are important 
differences in the spectrum. Two quasiparticle multiplets (particles and antiparticles) are 
four- fold degenerate, as in the SU(4) theory, but the six- fold degenerate multiplet of the 
SU(4) is split into a doublet with the mass M c _, a magnetic triplet with the mass M s t and 
a singlet mode with the mass M^ s . The ratios of the quasiparticle mass to the masses of the 
neutral modes are vastly different from the SU(4) Gross-Neveu model ratio M s /M e = \/2. 
They depend on the Coulomb interaction and the quasiparticle gap is much larger than the 
spin and the parity collective mode gaps. The gaps of the collective modes are not equal. 
In nanotubes where there are only two independent coupling constants we expect that 

2M C + M t + M s = (14) 
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(recall that since we deal with Majorana fermions, the masses may be negative, the spectral 
gaps being their absolute values). 



cos( v / vi : $i + - 1 ) 


= 0102) 


sin(v^F$( + )) 


= 


cos(v^F0( +) ) 


= Mi<7 2 , 


sin(v^Fei +) ) 


= 0"lA*2 


cos( v / vr$i _) ) 


= O-O03, 


sin( v / vr$i _) ) 


= ^0/^3 


cos(v^Fei- } ) 


= 


sin(V^©i _) ) 


= ^0^3 



IV. THE ORDER PARAMETERS OF DIFFERENT PHASES 

From now on we shall concentrate on the model describing carbon nanotubes. Differ- 
ent phases, shown in Fig. 3, correspond to different ground state phase lockings and are 
determined by the signs of the fermionic massesjsj], Eqs. pOjl . These sign changes are not 
reflected in the thermodynamics which is sensitive only to absolute values of the masses. 
The difference in the correlation functions may, however, be quite dramatic. 

We use the following conventions 



(15) 



(16) 

where {<7j} and {/Xj} (i = 1, 2, 3, 0) are order and disorder parameters of the 2D Ising models 
associated with the singlet (i = 1, 2, 3) and triplet (i = 0) Majorana fermions. A particular 
Ising model is ordered ((a) ^ 0) or disordered ((//) ^ 0) depending on the sign of the 
corresponding Majorana mass, M < or M > 0. In order to understand the structure 
and properties of the correlation functions, one has to recall that in the ordered phase of 
the Ising model, where (a) ^ 0, the correlation function {{n(u,q)n{—u,—q))) contains a 
coherent peak, while the correlation function of the ex's does not. 

There are six possible phases, two of them being Haldane spin liquids. Such liquids are 
characterized by the presence of a coherent triplet magnetic excitation (magnon) in the 
two-point correlation function of the staggered magnetizations. In non-Haldane disordered 
phases, the spectrum of triplet excitations is incoherent; however spin-singlet modes may 
be coherent. Table 2 and Fig. 3 show what operators acquire nonzero expectation values in 
the corresponding ground states. It will be assumed that in all phases the symmetric charge 
field is locked at <^ +) = 0, so that (cos v^4 +) ) °> ( sin v^ $ c +) ) = 0. 
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Table 2. 



Phase 


A 


(cos(a/vt$c~ 


} )) + (a a ) ^ (a ) 


B 


(cos(y / 7r$|f 


>)) ^ (/*«) 7^ ((To) ^ 


C 


(sin(-y/7r$c~ 


} )) ^ (/i a ) ^ K) ^ 


D 


(sin(y / 7T$c~ 


] )) ^ (ti a ) ± (/i ) ^ 


E 


(sin(A/7r$ c 


} )) ^ (<r a ) ^ (/io) ^ 


F 


(C0S( V /7T$C~ 


>)) ^ ^ (/i > ^ 




FIG. 4: The phase diagram. For each quadrant it is indicated which Ising models order and 
disorder parameters have non-zero ground state expectation values. The phases are separated by 
critical lines on which one of the particle masses vanishes. The green and red lines are Zi and 
SUiiQ critical lines respectively. The g\ = axis corresponds to the U(l) critical line. Phases A 
and D have density wave order. Under doping phases B, E and C,F develop a power law response 
to superconducting paring. 

Below we give a brief characterization of each phase; for more details see the subsequent 
discussion. 
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Phase A: g\ > \g^\. This phase has a long-range Bond Density Wave (BDW) order at 
T = (see Fig. 4 below). 

Phase B: < g\ < g%. This is a Haldane spin-liquid phase whose excitation spectrum 
contains a coherent triplet magnon with the mass M$ t t, associated with the correlation 
function of the site-diagonal spin operator S'~); see Eq. (jUJ). Under doping it develops 
a power-law response of the pairing susceptibility of type described further in the text. 
The corresponding superconducting order parameter transforms non-trivially under 
the lattice point group. 

Non-Haldane phase C: —g 3 < g\ < 0. The magnetic singlet mode with mass Ms, s 
becomes coherent in this phase. As a result a coherent peak appears in the correlation 
function of the Charge Density Wave order parameter A C dw, Eq. (JTHJ). Under doping 
this phase develops superconducting correlations in the s-channel. 

Phase D: g\ < — | £73 1 - There is a Charge Density Wave (CDW) long-range order in 
the ground state of this phase. Phases CDW and BDW are mutually dual. 

Phase E: g 3 < g\ < 0. This is another Haldane phase dual to phase B. Phase E has 
a coherent magnon in the site-off-diagonal spin density Eq. (j2f)j) . Under doping 

it develops superconducting correlations of the same type as phase B. 

Phase F: < gi < —g^. This phase is dual to phase C and has a coherent spin-singlet 
mode displayed by the BDW order-parameter correlation function. Under doping it 
develops superconducting correlations in the s-channel. 




o o o* 




FIG. 5: The dimerization in BDW phase 
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A. Order parameters, possible orderings and dominant correlations 



CDW. Already for geometrical reasons (the lattice is bipartite) one may expect that the 
system at half filling can develop a commensurate Charge Density Wave. Introducing locally 
averaged electron densities for a and b sublattices, 

pA r ) = ^ [i'li'uir) + iptip u (r + e!) + 4>l4> v (r - e 2 )] , v = a,b (17) 

we define the CDW order parameter (OP) as follows: 

A CDW (r) = p a (r) - p b (r) 

+ h.C. 

cr 

~ sin(v/^ +) ) cos(v^F$ c ~ ; ')<7i<72<7 3 (7o - cos( v / ^<I>c +) ) sin(v / 7r$ c _ ViA*2A*3A*o (18) 

This OP has a non-zero average value in phase D. In phase C where (sin( v /7r$ c ~' ) )) ^ and 
(fj, a ) 7^ (a =1,2,3), the most singular part of the CDW order parameter is proportional 
to /i . Since we are in the phase with (a ) ^ 0, operator p has a non-zero matrix element 
between the ground state and a state with one Majorana fermion. Therefore the correlation 
function of CDW OPs contains a coherent peak corresponding to an emission of the singlet 
magnetic mode with the mass Ms >s - 

BDW. The Bond Density Wave order parameter is similar to the CDW one, but is 
off-diagonal in the site indices: 

A BDW (r) = ^(r)[^(r + ei) - ~ e 2 )] - [V>l(r + ei) - ^(r - e 2 )]^(r) + h.c. 

a 

~ sin(v^F^ +) ) sin( v ^ : $JrViA*2A*3A*o + cos(^^ ] ) cos^^ - Vi^s^o- (19) 

This phase is dual to the CDW one. The BDW order parameter condenses in phase A. 
The dimerization ordering pattern in this phase is shown in Fig. 4. In the disordered phase 
F (which is dual to phase C), the spectral weight of the OP A BDW contains a coherent 
magnetic singlet mode with the mass Mg, s - 

Notice that the CDW and BDW OPs do not contain any oscillatory pieces. These show 
up in the density distributions that are not uniform across a sublattice but otherwise are 
consistent with the uniaxial symmetry of the nanotube: 

Pu(r) = \ fe^v(r + ei) + ^^(r - e 2 ) - ^„(r)l , v = a, b. 



2 
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Two incommensurate CDW OPs can then be constructed as follows: 
A 2w( r ) = Pa(r) + Pb(r) 



c 



(-) 



x sin( v / vr$i +) )yUi/i2Ai3° r o + icos( v / vr$ (+) )criO-2cr3/ i o + h.c; (20) 

A Ir7cw( r ) = P«(r) -p fe (r) 

„ e -*o*53 - LLL 2a ) + /i.e. 

~ K^e-^IR^ + L) X f] + h.c. (21) 

where k's are Klein factors. Since the dual antisymmetric charge field, is strongly 

disordered ((e*^ 9 * ) = 0), the OP A) 

ncow has zero expectation value in all sectors of the 
phase diagram. In phase F, where (a a ) ^ (a = 1,2,3), (/i ) ^ 0, the correlation function 
of the mCDW + OP displays a coherent peak corresponding to emission of a vortex particle 
with the mass M c _. 

The charge distribution corresponding to the OP A^^jy is depicted on Fig. El This 
kind of order can be induced by applying a modulated potential of sufficient strength that 
couples to A~ C£W . The potential must be strong enough to overcome the energy gaps of 
the corresponding excitation branches and to drive the system into a state with an induced 
order of the mCDW" type. 

Magnetization. The total magnetization is given by 

1 + 1 + 

= I R + 1 L + i{e 2 ^[( J R+ C rL 1 ) - {L+aR^} + h.c.} 
Its smooth part is equal to the sum of the chiral currents of the Majorana triplet: 

ja _ ja i to _ _ .i,al!>c 6 c 

R(L) — J 1,R{L) + J 2,i?(L) — — 2 Xr(L)Xr(L)- 



The oscillating part of the spin density is (Klein factors omitted) 



sin(V7r*c ))"oN + icos( v / ^<l > c >oN . (22) 



Here N = 0-203, 01/U203, 0102/U3) an d N = (01//203, ^102/^3, ^1/^2^3)- This OP is never 
coherent. 
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FIG. 6: The charge distribution for the mCDW order. 

On the other hand, to define the staggered magnetization S^ _ ^(r) associated with the 
bipartite lattice, we need to introduce locally averaged spin densities for the two sublattices 
(cf. Eq. (HZJ): 



s„ r 



1 r 

6 



iplaip v (r) + ipl(Tip v (r + ei) + ^cr^(r - e 2 ) 



v 



a, b 



(23) 



Then 



S<">(r) = i a (r)-s 6 (r) 

~ cos(v^F$[; +) ) cos( v / ^$[rVoN - sin(v^4 +) ) sin(v^F$ c ~ VoN 



R} 2 aL 2 ) + h.c. 



(24) 



The OP S^ _ ^(r) is coherent in the Haldane spin-liquid phase B, the corresponding particle 
representing a massive triplet magnon. Notice that in the expression ()24j) the vector field N 
)lays the role of the staggered magnetization of the effective antiferromagnetic spin-1 chain 



In full analogy with the BDW OP (|19|). one can construct the site-off-diagonal staggered 
magnetization: 

1 1 
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l - (i^IcrLx - R\<tL 2 ) + h.c. 



2 

~ cos v / ^*c +) sin \/7T$cTVoN + sin v^c^ cos V^$c~ )<7 oN (25) 

The structure of ()25|) indicates that phase represents a Haldane spin liquid and is dual to 
phase B. Indeed, in phase E it is the vector field N that can be regarded as the staggered 
magnetization of the effective S=l chain, and, hence, the spectral weight of the operator 
' contains in this phase a coherent triplet magnon. 

Pairing operators. Let us construct the s-wave superconducting (SC) order pa- 
rameter A s sc- To this end, we first build locally averaged, site-diagonal, singlet pairing 
operators: 

A ssc( r ) = ^ [Vv,TVv,i( r ) + VvrVVi( r + e + Vv,t^a( r - e 2 )] , v = a,b. 
The OP A s sc is then defined as 

A sSC (r) = A a sSC (v) + A b sSC (v) 

cr 

~ KiiK 2 ie~ l ^ e c ) (cos v^^c^i ^/^ + i sin v^^cT ViA^A^o) • (26) 

The amplitude of this OP acquires a finite value in non-Haldane phases C and F. These 
two phases, with coherent singlet magnetic modes in their spectrum at half-filling, exhibit 
power-law sSC correlations under doping. 

The s-wave type superconductivity is not the only possible singlet SC order one can 
imagine. One can introduce SC order parameter which transforms non-trivially under the 
point group: 

A m sc = j^at( r )ftM r ) -^H( r + e i) -VViO-e 2 )] - (T^l) 
~ ^2 a [L\ a R2- a + L2 a R\- a ) + oscillatory terms 

(7 

v ^ e ^ +> r c os( v ^$c~Vi/'2/'*30 r o + isin(v^F*crVi CT 20-3Aiol • (27) 



ro e v c 



This pairing OP exhibits power-law correlations in doped Haldane phases B and E. 

To see what kind of superconductivity this is, let us consider the lattice mean field 
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Hamiltonian: 






t(k) 


*•(*) 








-g*(k)A* 


g*(k)A* 






g(k)A 
-g(k)A 
t{k) 
t*{k) 



(28) 



where g(k) = — 1 + 2 cos(k x /2)e ly ^ ky / 2 . The spectrum is 

E 2 = \t(k)\ 2 +\g(k)\ 2 \A\ 2 . 
In the Haldane phase C under doping we get the following OP: 



(29) 



(30) 



The Aharonov-Bohm effect: the ^/-component of the vector potential is coupled to 
0R+ L Xa + Lf a R 1(7 ) - (RtL 2a + L+ a R 2(7 ) ~ 

cos( v / vr$c +) ) sin( v /7r$[r ) )o-i cr 2^30'o - smiy/ir^) cos(V^$c )^ij"2A*3/-«o (31) 

It is coherent in Haldane phase E with an emission of the magnetic singlet. 
The staggered flux around a hexagonal plaquette is 



$ = 6e i ^ 6 e 2iQx (L+L 1 - B+RJ + h.c. ~ {R fX f - L fX f )e 2iQx+ ^ /& + h.c. 



(32) 



where the flavor fermion is associated with $[» ^ field. Therefore there is an interesting 
possibility to drive a system into a peculiar critical state with a mixture of the magnetic 
singlet and nonmagnetic orbital mode by applying a magnetic field with the corresponding 
period. 



V. SMALL DOPING 



The large velocity difference between the symmetric charge modes and the other part of 
the spectrum holds a key to the stability of the approach at finite doping. For small doping 
kFd <C exp(— rcv/2e 2 ) the soliton mode of $g field, though becoming gapless, still lies above 
the others in the most of the momentum space (see Fig. 2b). This means that one can still 
integrate over and obtain Eq.® though with a prefactor < 1. Thus a small doping 
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will decrease the gaps. It may also give rise to finite decay rates for the collective modes. 
This is a pseudogap regime which is perhaps similar to the one existing in the underdoped 
state of the cuprate superconductors (see, for example, the recent paper The physics 

of this regime will be goverened by two energy scales: the scale of collective modes M c (see 
Eq.(|lU|) and the electronic scale M e . At energies smaller than M c the system is a Luttinger 
liquid, at M c < E < M e the collective modes contribute to all physical properties giving 
rise to strong enhancement of the magnetic susceptibility and specific heat. Above M e the 
effects of backscattering disappear. 

Let us briefly discuss a possible enhancement of the superconducting fluctuations in a 
doped regime. In those phases where the amplitude of the SC order parameter is frozen (see 
the discussion above) its power law correlations are determined by the scaling dimension of 
operator ex^p[i^/^^Q+ ]. This scaling dimension is equal to 1/AK C where K c is the renormal- 
ized Luttinger constant. At large doping the value of this constant is determined by the 
long-distant Coulomb interaction and is small, but close to the transition (small doping) it 
is close to one (see Q] for a more detailed analysis). Therefore there is a window of doping 
where Kq < 1/2 where the pairing susceptibility diverges. A more detailed description of 
the pseudogap regime will be given elsewhere. 



VI. CONCLUSION 



Let us summarize our results. We have analyzed the effect of the long-range Coulomb 
interaction on the low-energy properties of the armchain nanotube. We have shown that 
the spectrum of the system is very rich. The phase diagram includes six phases differing 
from each other by their response functions. Under doping four of them have enchanced 
superconducting fluctuations. 

SWCNs exhibit two types of behavior: either band insulators or Luttinger liquids. Arm- 
chair CNs are believed to be metallic and be of the Tomonaga-Luttinger (TLL) type. This, 
in fact, is only true away from 1/2 filling. As shown in this paper, the picture changes dra- 
matically at 1/2 filling when effects of strong correlations become crucial. To observe these 
effects in experimental conditions, the chemical potential must be fine tuned to the value 
/i = at which the Fermi "surface" is represented by two degeneracy (Dirac) points. Then, 
due to the special role of Umklapp processes in the presence of unscreen Coulomb interac- 
ts 



tion, on lowering the temperature, one must be able to see a crossover from the metallic, 
TLL behavior to an insulator discussed in the paper. Such insulating state is Mott-type, 
with a much richer spectrum than in the Hubbard model. One important difference is that, 
unlike in the Hubbard model, the spin excitations have spectral gaps. 

The estimates of interaction matrix elements for carbon nanotubes are provided in jjj: 

e 2 

9i = ^We i = c ' t ' s 

ft ^0.4, ft « 0.5, ft « -1.3 (33) 

where N is the number of transverse bands of the nanotube and e m 1.4 is the dielectric 
constant. This estimates indicate that realistic nanotubes belong to phase B. Eq. fjlOj) also 
gives the following estimates for the spectral gaps of the collective modes: 

im -i = ^iLl < 34) 

Since a is the short distance cut-off of the bosonic theory and is therefore non universal, 
this formula contains a certain degree of ambiguity. For noninteracting electrons ira = a 
(the lattice spacing). Taking N = 10, v/v c = 1/3 and a = 0.246nm we get the estimate of 
order of 0.1 ev. The single particle gaps should be even greater. 

Since the gaps should be quite sizeable, one may wonder whether they have not been 
already observed. The existing techniques produce simultaneously carbon nanotubes of dif- 
ferent sizes and chiralities and one has to select the relevant ones using some criteria. The 
simplest one is to differ between metallic and insulating (semiconductor) nanotubes. Natu- 
rally, on the first glance Mott and band insulators look alike. Therefore it is possible that 
the Mott insulating armchair nanotubes have been overlooked being taken for semiconduct- 
ing ones. To distinguish Mott insulator from a band one one has to measure transport and 
magnetic properties and compare the gaps. The clearest sign of 'mottness' is the difference 
in gap sizes in different response functions. This is the feature to look for experimentally. 
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APPENDIX A 



The original Hamiltonian describing noninteracting electrons on a honeycomb lattice can 
be compactly written in the two-sublattice (Nambu) representation: 



H = ^Mk)^, 



(Al) 



where 



V>a(k) 



W(k) 



V 



t(k) = 1 + 2 cos(A; :E /2)e i(v ^ /2) ^ 




(A2) 
(A3) 



and the sum in (jAlj) goes over the Brillouin zone. The spectrum has two degeneracy points 
(nodes) at Qi^ = (±47r/3, 0). Linearizing the noninteracting Hamiltonian near these points 
yields two cones associated with (2+l)-dimensional massless Dirac fermions: 



ft(Qi,2 + p) = v(r y p y =F r x p x 



V3t/2. 



(A4) 



When a two-dimensional sheet of graphite is wrapped to produce an armchair nanotube, 
k y gets quantized, and the lowest-energy subband correspoding to k y = stays gapless. 
The resulting problem is one-dimensional because the wave function does not depend on y. 
Projecting the fermionic annihilation operators of the a and b sublattices onto the k y = 
subspace, we get: 



(A5) 











( r 2 (x) \ 


v A(r) ) 




{ h(x) ^ 







The effective ID Hamiltonian is brought to its canonical diagonal form 

H Q = -iv / d ^ (R)d x Rj - L]d x Lj 



(A6) 
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by linear transformations: 



n = (Rx + LJ/V2, l 1 = {-R 1 +L 1 )/y/2 
r 2 = {R 2 - L 2 )/V2, h = {R 2 + L 2 )/V2 



(A7) 



The structure of these transformations reflects the fact that different Dirac points, k x = Q 
and k x = —Q, are not associated with fermions of different chiralities, as it is the case for 
standard chains and ladders; instead each of these points is characterized by a pair of right 
(Rj) and left (Lj) fields. This is because the gapless spectrum of the armchair nanotube 
keeps the memory of the two-cone Dirac structure of the dispersion law in the 2D graphite 
(in fact, the ID spectrum is obtained from the 2D one by cutting the two cones by the plane 
k y — 0). For this reason, in case of nanotubes, smooth components of the physical fields are 
contributed not only by the diagonal "currents", R^Rj and LjLj, but also by off-diagonal 
"mass bilinears", RjLj and LjRj. 

The chiral fermionic fields can be bosonized in terms of chiral bosonic fields $^ L : 



Here a is the short-distance cutoff in the bosonic theory, and Klein factors obeying 

the algebra {^j a ,^j'a'} = Sjj/8 aa /. The product of the four Klein factors, T = k^ki^k 2 ^k 2 i, 
satisfies T 2 = 1. Since T is not a dynamical variable, we can conveniently choose T = 1. 

In the bulk of this paper, we adopted the description in terms of four scalar fields, 
$^, ± ' ) , and their dual counterparts, 8^, known from earlier studies of the two- 

channel Kondo problem . These fields describe the symmetric and antisymmetric charge 
excitations (equivalently, the "charge" and "flavor" modes), 



as well as the symmetric and antisymmetric spin excitations (or the "spin" and "spin-flavor" 
modes), 




(A8) 



$ i ±} = 2 ( $ 1T + $ U ± $ 2 T ± $21) 



(A9) 




Here $ JO - = + The corresponding dual fields, and 0^ are obtained from the 
above expressions by replacing <3>j CT by 6 Jcr = —&f a + &j a - 
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APPENDIX B 



The interaction can be written as 

[p a + Pb]{ri)V aa (r 12 ){pa + Pb](r 2 ) + 

2p a (r 1 )[V ab (r 12 ) - V aa {r 12 ))p b (r 2 ), (Bl) 

where, in the low-energy limit, the local electron densities on the a and b sublattices are 
represented by 

p a (r) -> p r {x) + e 2iQx M r (x) + e~ 2iQx M+(x), 
p b {r) -> Pl {x) + e 2iQx Mi{x) + e~ 2iQx M+{x), 



p r = rjrx + r+r 2 , pi = Ijh + ql 



2 



M r = rtn, M l = lth (B2) 
Substituting Eq. (jB2|) into Eq. (|Bl|) and dropping the oscillatory terms we get four terms: 

• (p r +pl)lV aa (r 12 )(p r + pi) 2 , 

• 2[U ab (0) - U aa (0)}p r (x) Pl (x), 

• 2U aa (2Q)[M r (x)M+(x) + M ; (x)M+(x)], 

• 2U ab (2Q)(M r Mf + M t M+), 

where U(0) and U (2Q) stand for the Fourier transforms of the interaction potentials. The 
first term here gives the Luttinger coefficient renormalization. The 'backscattering' inter- 
action expressed in terms of the standard Dirac fermions with flavour indices 1,2 looks as 
follows: 

- gi (R+ a L 1>cr - RlMv + h.c.) 2 + 
g 2 [(L+ a L 1 , a )(R+ y R 2 , a ,) + h.c.] + 

9z{L 2i0; Rl,a - R 2 ,aLl,a){LiyR 2 y - RtyL 2 y) (B3) 

where the couplings are expressed in terms of the Fourier components of the interaction: 

91 = U aa (0) - U ab (0), 

92 = U ab (2Q)-U aa (2Q), g 3 = U aa (2Q) + U ab (2Q). (B4) 
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In the UV model, which apart from the Hubbard (on-site) interaction U also includes the 
interaction between electrons on nearest-neighbor sites, V, the couplings gi are given by 



9l = U-3V, g 2 = -U, g 3 = U 



(B5) 



The bare Hamiltonian for R, L is standard: 




(B6) 



As already explained in section I and II, in the presence of long-range Coulomb potential, 
Umklapp processes with the structure R\R\L^Li + h.c. represent the strongly relevant part of 
interaction (jB3|) . In the bosonic language, these are the processes containing the symmetric 
charge field Using bosonization rules (jA8|) - (jAlOj) . one straightforwardly derives Eq. 
(0). The g 2 -part of (|B3|) does not contribute. 

If the interaction is the Coulomb one, then the couplings #1,(73 are positive. It also looks 
likely that g 3 > g\ and \g 3 — gi\ « gi j3 . At these circumstances we have to consider two 
possibilities: (i) g\ < g% (Haldane SL) and (ii) g\ > g% (CDW) . 



[1] R. Egger and A. O. Gogolin, Phys. Rev. Lett. 79, 5082 (1997). 
[2] C. Kane, L. Balentz and M. P. A. Fisher, Phys. Rev. Lett. 79, 5086 (1997). 
[3] H. Yoshioka and A. A. Odintsov, Phys. Rev. Lett. 82, 374 (1999); Phys. Rev. B 59, R10457 
(1999) . 

[4] L. Levitov and A. M. Tsvelik, Phys. Rev. Lett. 90, 016401 (2003). 

[5] L. Glazman et al, Phys. Rev. B45, 8454 (1992); H. J. Schulz, Phys. Rev. Lett. 71, 1864 



[6] R. Saito, G. Dresselhaus, and M.S. Dresselhaus, Physical Properties of Carbon Nanotubes, 

Imperial College Press, 1998. 
[7] C. Wu, W. V. Liu and E. Fradkin, Phys. Rev. B68, 115104 (2003). 

[8] D. Shelton, A. A. Nersesyan and A. M. Tsvelik, Phys. Rev. B 53, 8521 (1996) (1996); A. O. 

Gogolin, A. A. Nersesyan and A. M. Tsvelik, Bosonization in Strongly Correlated Systems, 

Cambridge University Press (1999). 
[9] D. G. Shelton and A. M. Tsvelik, Phys. Rev. B53, 14 036 (1996). 



(1993). 



23 



[10] I. V. Krive, A. A. Nersesyan, M. Jonson and R. I. Shekhter, Phys. Rev. B52, 10865 (1995). 
[11] H. H. Lin, L. Balentz and M. P. A. Fisher, Phys. Rev. B58, 1794 (1998). 
[12] Al. B. Zamolodchikov, Int. J. Mod. Phys.AlO, 1125 (1995). 

[13] U. Schollwock, S. Chakravarty, J. O. Fjaerestad, J. B. Marston and M. Troyer, Phys. Rev. 

Lett.90, 186401 (2003). 
[14] E. Papa and A. M. Tsvelik, Phys. Rev. B63, 5109 (2001). 
[15] V.J. Emery and S. Kivelson, Phys. Re, B 47, 10812 (1992). 



24 



